Inertial particles driven by a telegraph noise 
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We present a model for the Lagrangian dynamics of inertial particles in a compressible flow, 
where fluid velocity gradients are modelled by a telegraph noise. The model allows for an analytic 
investigation of the role of time correlation of the flow in the aggregation-disorder transition of 
inertial particle. The dependence on Stokes and Kubo numbers of the Lyapunov exponent of particle 
trajectories reveals the presence of a region in parameter space (St, Ku) where the leading Lyapunov 
exponent changes sign, thus signaling the transition. The asymptotics of short and long-correlated 
flows are discussed, as well as the fluid-tracer limit. 

PACS numbers: 47.52.+j, 47.55. Kf,47.27.eb 



INTRODUCTION 

The spontaneous formation of clusters of particles sus- 
pended in chaotic flows may originate from two different 
physical processes: compressibility of the fluid flow and 
particle inertia. In the first case particles are trapped 
in regions of ongoing compression, while in the second 
case inertia causes their ejection from vortical regions. 
The underlying link between these two phenomena is 
manifested in the limit of weak inertia, in which par- 
ticle dynamics in incompressible flows can be approxi- 
mated bythat of tracers in a weakly compressible veloc- 
ity field [IB- 

Clustering processes have been extensively studied, 
both in the case of compressible flows [!, 0, H, [|| and iner- 
tial particles @, 0, H, [H • Fractal dimension of the clusters, 
obtained from the ratios of Lyapunov exponents of parti- 
cle trajectories, has revealed a powerful tool to quantify 
the intensity of the clustering 0, 0] , and has motivated 
further studies of Lyapunov exponents of inertial parti- 
cles 

Aggregation is an extreme form of clustering, which 
occurs when the trajectories of different particles tend to 
point-like clusters. In this situation the senior Lyapunov 
exponent of particle trajectories become negative, signal- 
ing that nearby particles do not separate exponentially 
as expected in a chaotic flows, but their paths eventually 
coalesce. This phenomenon has been described for fluid 
tracers in compressible flows [HI, EH, as well as for light 
particles in incompressible flows [7| . 

Surprisingly, it has been shown that in compressible 
flows, where fluid trajectories coalesce, large enough in- 
ertia can induce a transition from the strong clustering 
regime into a week clustering one, where particle tra- 
jectories remains chaotic and the senior Lyapunov ex- 
ponent is positive [l3|, 14 , [l5|, [l6| . Analytic results on 
this aggregation-disorder transition have been obtained 
under the assumption that the carrier flow is short- 
correlated in time. 

It is therefore natural to ask how time correlations of 
real turbulent flows can influence this phenomena. As 



pointed out by recent numerical studies of fluid tracers 
in compressible flows [fjj], time correlations are responsi- 
ble for an increase in the level of compressibility required 
for the transition to the strong-clustering regime. In the 
case of inertial particles in has been found [lCJ that time 
correlations cause increase the (positive) Lyapunov ex- 
ponent (i.e. make particle flow more chaotic) already at 
weak inertia, which can not be pred icted in the frame- 
work of short correlated flows [9l.rT7, 3, 19j |. 

Here we discuss a simple time-correlated model flow, 
recently introduced in [201 ] . which allows to get a deeper 
insight on the interplay between inertia and compress- 
ibility, and to investigate analytically the dependence on 
Stokes and Kubo numbers of Lyapunov exponent of par- 
ticle trajectories. 



MODEL 

The dynamics of two small inertial particle, whose den- 
sity is much larger than the fluid density is dominated by 
the viscous drag. Hence, the equation for their separation 
R(t) = Xi(<) - X 2 (t) and relative velocity V reads 



R = V , V 



1 



(V-AU(X 1 ,X 2 ,t)) (1) 



where r is the Stokes time of the particles, and AU is 
the difference of fluid velocities at particle positions. 

To investigate the behavior of Lyapunov exponent one 
has to consider separations smaller than the viscous 
length-scale of the flow. At these scales the fluid ve- 
locity difference can be written in terms of the velocity 
gradients Sij = djUi as AUi = SijRj. The Lyapunov ex- 
ponent A is determined by the contributions of the local 
gradients Sij(t) experienced by the two nearby trajec- 
tories. Note that the statistics of velocity gradients in 
the reference frame of an inertial particle differs from the 
Eulerian statistics due to preferential concentration. 

Let us now consider an idealized one-dimensional flow, 
in which the fluid velocity gradient in the reference frame 
of a particle is modelled by a telegraph noise s(t), i.e. a 



2 



noise that switches randomly between two fixed values 
±s. Stokes number can be defined in term of the fluid 
gradient intensity as St = st. Equations |T]) reduce to 



v, v = -^ + s M R 

T T 



(2) 



When s(t) = +s the system is expanding, and evolves 
towards the asymptote V = aR, where 



1 



a= — (-1 + VT 
zt 



4st) 



(3) 



The particles separate exponentially with smaller rate 
than fluid trajectories a < s. When s(t) = — s the system 
is contracting, and two different scenarios are observed, 
according to the relative intensity of fluid gradients and 
inertial drag. For small Stokes numbers St < 1/4 the 
system evolves toward the asymptote V = PR, where 



1 



P = — (-1 + VI - 4st) for st < 1/4 . 
2r 



(4) 



Particles slow down less efficiently than the fluid, and 
hence their separation goes to zero with a faster expo- 
nential rate \j3\ > s. For large Stokes number St > 1/4 
the system has two complex conjugate eigenvalues 



P 



(-1 ± iy/AsT- 



1) for st > 1/4 



(5) 



and the solution decays exponentially with a clockwise 
spiral motion in phase space (-R, V). This means that 
particles can cross the R = axis, i.e. collide, with 
non-zero relative velocity, giving orig in to a shock i.e. 
discontinuity in the velocity field [22j. Notice that the 
solution is properly defined also for R < 0. The change 
of the sign of R can be interpreted as the fast particle 
overcoming the slow one. 

Shocks start to appear for the critical value St* = 
1/4. The presence of such critical Stokes number is 
characteristic of flows where fluid velocity gradient are 
bounded, and its value is determined by the intensity of 
the strongest negative gradient (here — s). Conversely, 
if the statistics of fluid velocity gradient is unbounded, 
shocks can appear for arbitrarily small Stokes number, 
but the y ar e exponentially suppressed in the limit St — > 

@M 11. 

The transition rates v\ from — s to s and v-i from s 
to — s determine the fraction of time spent by the parti- 
cle pair in regions of ongoing compression and expansion 
(respectively, vij(y\ + v<x) and v\j(v\ + i^))- The mean 
value of the noise s(t) is 



s = (s(t)) = -sAv/v, 



(0) 



where v = v\ + v-i and Av = V2 — v\. Noise fluctuations 
s(t) = s(t) — sq are exponentially correlated: 



(S(t)S(t')> 



As 2 V\V2 



exp(-^|< - t'\) 



(7) 



The ratio between the correlation time of the flow \Jv 
and the Lagrangian separation time 1/s defines the Kubo 
number Ku = sv^ 1 . 

Notice that in potential flows both fluid tracers and 
inertial particles spend more time in regions of a local 
compression than in expanding regions, and hence Av > 
0. The actual value of Av can be determined by assuming 
statistical homogeneity and isotropy of the flow. Thanks 
to these symmetries, the mean position of a particle is 
solely determined by its initial position (X(t)) = Xq and 
hence the distance R is statistically conserved. 

The peculiar nature of the telegraph noise, namely 
the fact that its square is deterministic, allows to ob- 
tain closed equations for (R(t)). Averaging over different 
realizations of the noise, and using the derivation for- 



mula [24 1 



dt 



<6 s *itt)> 



(8) 



which hold for a generic noise with zero mean and cor- 
relation (£t£e) ~ exp(— v\t — t'\), one obtain the linear 
system 



d 1 
y = -My 



dt 



(9) 



where 
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~{1 + vt) 



(10) 



This linear system allows for a constant solution (R) = 
R when sqv(1 + ut) + s 2 = which together with Eq. (JB]) 
gives 



Av = 



1 



(11) 



Notice that inertia reduces the trapping of particle in 
compressing regions. This effect is encoded in the behav- 
ior of the parameter Av, which is a decreasing function 
of Stokes time t, and vanishes in the limit r — -> oo. 

The positivity of v\ requires Av < v, and gives an up- 
per bound for the Kubo number achievable in the system: 



Ku < Ku* 



2St 



VTTTSi-1 



(12) 



For longer correlation time particle are irreversibly cap- 
tured by contracting regions. In the limit of fluid tracers 
r — > the constraint becomes Ku < 1. The accessible 
region in (St,Ku) parameter space is shown in Figured] 



3 




0.01 



0.01 



0.1 



1 

St 



10 



100 




FIG. 1: Parameter space (St,Ku). Shape transitions in the 
probability density function (pdf) of x occur for St = St* = 
1/4 (vertical solid line), rh — 1 (dashed line) and m = 1 (dash- 
dotted line). Labels refers to the pdfs shown in Fig. [2] [3] Grey 
region is not accessible by the model because Ku number does 
not fulfill condition (|12[) . Crosses represent the boundary of 
the chaotic region. 



STATISTICS OF PARTICLE- VELOCITY 
GRADIENTS 



FIG. 2: Pdf of x in the small Stokes regime (St = 1/8) for 
different value of Kubo number: K — l/8(a), Ku = 1/16(6) 
and Ku = l/24(c). 

\x\ gives the probability of strong particle velocity gradi- 
ents, and is therefore related to the probability of shocks. 
Different solutions are found in the small-Stokes-numbcr 
and large-Stokes-number regime. 



From the system ^ one can obtain a closed equation 
for the particle- velocity gradient a = V/R: 



Small Stokes number 



& = -a 2 - T^ 1 [a - s(t)] . 



(13) 



The Lyapunov exponent of particle trajectories can be 
obtained from Eq. (|13p as an ensemble average (a) over 
different realizations of the noise. Equations for the sta- 
tionary distribution of x = a+l/ (2r), are easily obtained 
with the same procedure adopted to derive Eq. (|9ll0p : 



4r 2 



so 

7' 



i 

17 



so , 

7' 



q x = 

(14) 
vq = 

(15) 



where q(x) = J sp(x, s)ds/r and p(x, s) is the joint sta- 
tionary pdf of x and s. From the first equation one gets 
q = — (C + (l/4r 2 — x 2 + so/r)p) where C is the mean 
flux of x, and finally 



4t 2 



( 4 — ) (4^ 



2 \ S 

x — v — 



(16) 
C(2x -v) = 



The asymptotic behavior p ~ Cjx 2 (C > 0) for large 



When St < St* the unique positive integrable solution 
of Eq. (fT6|) is obtained under condition of zero flux (C = 
0): 



p=C_ 



(w — x) m 1 (x — w) T 



(w + x) m+1 (x + w) rh+1 
and zero otherwise, where 



x G (w, w) (17) 



w = vl + 4st/2t m 
w 



sf\ - Ast/2t m 



(v + Av)/4w 
(v - Av)/4w, 



(18) 



The solution is localized in the compact interval (w, w). 
Its shape is determined by the values of m and rh. For 
m < 1 (low frequency) the pdf is peaked around the two 
border values. For m < 1 < rh (intermediate frequency) 
it vanishes at w, and finally when m > 1 (high frequency) 
it vanishes both at w and w (See Fig. [3]). Indeed, when 
St < St*, the solution of the linear system for (R,V) 
oscillates between two asymptotes V = (w — \/2t)R and 
V = (w — 1/2t)R according to the sign of the noise. If the 
noise frequency is low the system has enough time to get 
close to the two asymptotes and the pdf is peaked around 
them. Conversely, when the sign of the noise switches 
frequently, the system does not have enough time to reach 
the asymptotes and oscillates rapidly around the mean 
value. 
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FIG. 3: Pdf of x in the large Stokes regime (St = 1) for for 
different values of Kubo number: Ku = 0.1(d), Ku = 1(e) 



Large Stokes number 

When St > St* the solution of Eq. (fTTH) consists of two 
different parts: 

p{x) = C Pl (x) + G 2 P2{x) (19) 
The first one is the solution of Eq. (p~6|) with C = 1, 



Pi 



^n(arctan(:E /w)) 

\w + x\ m+1 (x 2 + w 2 ) 

\w + y\ m+1 {v-2y) e _ w(arctaa(y/a)) 
\y — w\ m ^ 1 (w 2 — y 2 ) 



(20) 



while the second one is the right tail of the solution of 
the homogeneous (fluxless) equation: 



P2 



\w + x\ m+1 (x 2 +w 2 ) 

n(arctan(a;/ij})) j I \ 



(21) 



Here I Wi00 (x) is the characteristic function (indicator) 
of the interval (w, oo), and 



w = yisr +T/2t to = (i/ + Ai/)/4to 
il) = V4sr - 1/2t n = (v — Av)/2w. 

Conditions for determining C and C2 are 

p(x)dx = 1, / q(x)dx = . 



(22) 



(23) 



These conditions guarantee the continuity of the pdf in 
the limit St — * St* (see Appendix ) The pdf obtained 
in the regime St > St* is extended over all real x, 
with power-law tails p ~ C/x 2 for large |x| that are 
due to shocks, which occur for large negative values of 



s(t). For short-correlated noise (to > 1), the pdf is 
characterized by an asymmetric core localized between 
— w and w. When to < 1 a singular peak arises at 
x = w (See Fig. [3]). This behavior is easily understood in 
term of solution of the linear system for (R, V). When 
s(t) = +s the solution converges toward the asymptote 
V = aR = (w — 1/2t)R. This produces the peak at 
x = w, provided that the correlation time of the noise 
is long enough (to < 1) to get close to the asymptote. 
The large-Stokes number regime is hence characterized 
by an infinite series of shocks alternated to 'quiet' phases 
in which the particle- velocity gradient relaxes toward a. 

The schematic in Figure [T] summarizes the different 
regions in the parameter space (St, Ku) where shape 
transitions occurs in the pdf of x. 



LYAPUNOV EXPONENTS 

The Lyapunov exponent of incrtial-particle trajectories 
can be written in terms of the mean value x as A — 
x — 1/2t. For St < St* the mean value x can be written 
as 

(w — w)rh 

x = w 



(to + to) 

F\ (to + 1, m + 1,to+ 1,to + to + 1, u, v) 
Fi(rh+, m + 1, fh + 1, rh + m, u, v) 



(24) 



where 



w — w 
w + w' 



w 



2w 



(25) 



and F\ is the hypergeometric function of two 
variables [2^ . The behavior of the Lyapunov exponent as 
a function of Stokes and Kubo number is shown in Fig- 
ure |U Lyapunov exponent is negative for small Stokes, 
and decreases approximatively as As -1 ~ —Ku at in- 
creasing Ku numbers. In the limit St —> it recovers 
the actual value for fluid tracers 12011: 



Aq = lim A = —s 2 jv 



(26) 



Notice that fluid tracers are always in the aggregation 
regime within this model, as signaled by the negative 
value of Aq. A sharp negative minimum A = —2s is found 
for St = 1/4 at Ku = (V2 + l)/2. It corresponds to 
the maximum aggregation of particle. As the correlation 
time of the flow decreases, the minimum becomes less 
pronounced, and it moves to larger Stokes numbers. A 
region of positive Lyapunov exponents, is present in the 
parameter-space for St > 4.3. The isoline of vanishing 
Lyapunov exponent which border this region represent 
the transition from the strong clustering regime (A < 0) 
to the chaotic regime (A > 0). Indeed, Figure [5] shows 
that, as Stokes number increases, the interval of Kubo 
numbers appears where the Lyapunov exponent grows, 
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FIG. 4: The isolines of the Lyapunov exponent A on the plane 
of Stokes and Kubo numbers. The isolines are spaced every 
0.02s The boundary of the chaotic region (A > 0) is repre- 
sented by the dashed line. 



FIG. 5: The Lyapunov exponent A as a function of the Kubo 
number. An interval of positive Lyapunov exponents appears 
for St > 4.3. The dashed line represents the asymptotic be- 
havior for Ku — Ku* . 



and eventually becomes positive. This can be understood 
as follows: to achieve an effective mixing the correlation 
time of the fluid gradients must be long enough to pro- 
vide substantial stretching of particle trajectories, but 
not too long, to avoid particle segregation in compress- 
ing regions. Therefore the chaotic region is confined in 
a window of Ku numbers between a lower and an upper 
bound determined by stretching efficiency and particle 
trapping respectively. 

Let us now consider the behavior of the Lyapunov mo- 
ments 7„, defined as (R n ) ~ exp(7„t). The evolution of 
(R n ), for n positive integer, is determined by a closed 
system of 2(n + 1) equation: 



A and B are tri-diagonal matrices with the elements 



-{R n - l V l 
dV 



(n - i)(R n - l - 1 V l+1 ) 
i 



T 

i 



rti' 



sR n ~ l V l ) = 



(R n -' l V l ) 
-(5 + s )vH n ~ <+1 V* -1 ) 

T 

{n-i){sR n - i - l V i+1 ) 



- + 

-soiSR^^V 1 - 1 ) 

T 



(27) 



(28) 



where i = 0, 1, • ■ ■ , n. The n — th Lyapunov moment j n 
is the largest solution of the 2(n + l)-th equation: 



det 



A C 
D B 



= 0, 



(29) 





7n + i/t 


i = 1 , n 






= -iso/r 


i = 2,n 






= i — n 


i = 1 , n — 


1 


Bis = 


7„ + v + i/t 


i = 1 , n 




Bij-i 


= isq/t 


i = 2,n 






= i — n 


i = 1 , n — 


1 



(30) 



and C and D are sub-diagonal matrices with the ele- 
ments CVi-i = —i/r for i = 2,n and Di.i-i = i(s§ — 
s 2 )/r for i = 2,71. In the limit St — > one recovers the 
Lyapunov moments of fluid tracers [20| : 



7 " = \/© 2 + s2(n2 " n) "i (31) 

Notice that for fluid tracers the sign of the separation R is 
preserved, and hence the exponents 7„ coincide with the 
exponents 7„ defined by (|-R| n ) ~ exp("f n t). The asymp- 
totic linear behavior of "f n for large n is the hallmark of 
the presence of an upper bound for velocity gradients. 



Short-correlated flows 

Let us now examine the limit of short-correlated flow. 
The limit v — > oo must be taken keeping constant 
s 2 /v = D, in order to recover (^-correlated noise fluc- 
tuations (s(t)s(t')) = 2D6(t — t') In other words, while 
Kubo number Ku = sjv tends to zero, Stokes number 
St = st must grow, so that the product 



KuSt 



l jvT = DT 



(32) 
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FIG. 6: Lyapunov exponent A in the short-correlated limit. 



remains constant. In this sense the short correlated 
limit for inertial particles correspond always to the large- 
inertia case. In the delta-correlated limit the relevant 
time-scale associated to fluid velocity gradients is given 



by the Lyapunov exponent of fluid tracers Aq 



7" 



—D. This is confirmed by the collapse of particle Lya- 
punov exponents in the short-correlated limit once their 
intensity and Stokes times are re-scaled with |Ao| (see 
Fig. [6J. A noticeable minimum is observed for |Ao|r = 
KuSt ~ 0.05 and a transition to chaos, i.e. from nega- 
tive to positive A, occurs for KuSt > 1.6. These fea- 
tures are in qualitative agreement with previous ana- 
lytic and theoretical results obtained in the framework 
of 5-correlated flows 14, 18| ■ Notice that in those studies 
Gaussian statistics is assumed for velocity gradients, at 
variance with our model in which only the two values ±s 
are allowed. Therefore quantitative details such as the 
exact position of the minimum can be different. 

We remark that in the short-correlated asymptotics, 
the fluid-tracers limit becomes singular. Indeed one has 



lim so = lim — 



u{l + vt) 



-D r = 
r > 



(33) 



which signals that fluid tracers are still preferentially 
attracted by regions of ongoing compression, while in- 
ertial particles with arbitrary finite r are not. Notice 
that in the limit v — ► oo fluid velocity gradients be- 
comes unbound because s — ► oo, and for fluid tracers 
we recover the quadratic behavior of Lyapunov moments 
7„ = D(n 2 -ri). 



Long-correlated flows 

Time correlation of fluid gradient is bounded by the 
condition (fT2")) . When Ku = Ku* the transition rate v\ 



FIG. 7: Lyapunov exponent A for Ku = 1 (solid line) and 
on the line Ku = Ku* (dashed line), corresponding to the 
longest correlation of velocity gradients achievable within the 
model. Upper inset: asymptotic behavior for large Stokes 
number As^ 1 ~ St" 2 '^ 3 ', {Ku = 1). Lower inset: asymptotic 
behavior for small Stokes number A — Aols -1 ~ St'' . 



from 
s to 



-s to s vanishes, while the transition rate from 
s reaches its maximum v* 



(Vl + 4sr-l)/(2r). 
Particles initially seeded in expanding regions are grad- 
ually captured by contracting ones, where they remain 
trapped forever. Therefore the population of expanding 
regions decreases exponentially as 



P(t) ~ exp(-i/*f). 



(34) 



Moments of particle separations will hence evolve asymp- 
totically according to 



\R\ 



e nat P(t) + (1 - P{t))e nRe(l3)t f{Im{P)t) (35) 



where (3 is given by ((5]) and f(t) is 2tt- periodic function. 
Notice that v* = a and therefore for n = 1 the decreas- 
ing fraction of particle in expanding regions is exactly 
balanced by the exponential growth of their separation. 
From (|34l35p one gets the Lyapunov moments 



7n 
In 



a(n-l) for n>{a- Reifi))- 1 
-- Re(P)n for n < (a - i?e(/3))" 1 



In the limit of fluid tracers one gets 



7„ = s(n - 1) for n > 1/2 
7„ = — sn for n < 1/2 



(36) 



(37) 



in agreement with Eq. (|3T|) . 

Finally, the Lyapunov exponent along the critical line 
Ku = Ku* is A = {d%/d)\ n=0 = Re{/3) In Figure we 
compare its behavior with that of the Lyapunov exponent 
along the line Ku = 1. For small Stokes number both 
of them recover the fluid-tracers limit A — > Aq = — s but 



7 



with different power law behavior. On the line Ku = 
Ku* we have |A — Aols" 1 ~ St, while for Ku = 1 we have 
|A — A |s _1 ~ St 2 (see lower inset of Fig. [7]). For large 
Stokes number Lyapunov vanishes as As" 1 ~ —St -1 on 
the line Ku = Ku* and as As" 1 ~ St~ 2 / 3 for Ku = 1 (see 
upper inset of Fig.|7J). In between these two asymptotics 
a sharp minimum appears for St = St*. 

Notice that the asymptotic decay St~ 2 ^ 3 here shown 
for long-correlated flows, have been already predicted 
and observed also for 5-correlated flows [17|, |l9j. The 
agreement between these results confirms that in the 
large-Stokes number asymptotics the role of time corre- 
lation becomes negligible and particles behave as if sus- 
pended in (S-correlated-in-time flows. 



CONCLUSIONS 

We discussed the Lagrangian dynamics of inertial par- 
ticles in a simple time-correlated compressible flow, in 
which fluid velocity gradients in the reference frame of 
the particle are modelled by a one-dimensional telegraph 
noise. In spite of its simplicity, the model allows to 
take into account consistently the different Lagrangian 
weights of region of ongoing compression and expansion, 
and it reproduces the phenomenon of trapping of parti- 
cles in compressing regions for long-correlated flows. 

The peculiar nature of the telegraph noise allows to 
investigate analytically the effects of time-correlation of 
velocity gradients on the chaoticity of particle trajecto- 
ries, and to study the dependence on Stokes and Kubo 
numbers of Lyapunov exponents. We discussed both the 
asymptotics of long and short-correlated flows, as well as 
the fluid-tracers limit. 

For large Stokes number, a regime characterized by the 
formation of shocks, we found a chaotic region in parame- 
ter space (St,Ku), where the leading Lyapunov exponent 
becomes positive. Inertia is therefore responsible for a 
transition from a strong clustering regime, originated by 
the compressible nature of the flow, to a chaotic regime. 
The latter is observed in a range of Ku numbers such that 
the time correlation of fluid gradients is long enough to 
provide substantial stretching, but not too long to get 
particles trapped in compressing regions. 



APPENDIX A 



To prove the continuity of the pdf for St = St* we first 
notice that the limit St — > St* is equivalent to w —* 0. 

The limit w — ► of the solution (fT7)l in the small- 
Stokes number regime is p[x) = Cif(x)Io, w {x), where 



/(*) 



(w — x) r 



(w + x) m+1 x 2 



-{v-Av)/2x 



(38) 



and Iq. w (x) is the characteristic function of the interval 
(0,w)' 

To study the limit of the solution (fl9|) let us consider 
the limits ofpi(x) and P2(x) when w — > 0. The latter is 
easily obtained as: 



p 2 (x) ~ e^-^^fix) 



(39) 



Now let us rewrite p\ {x) as 
\w-x\ m ~ Y 

Pi = 



-n(arctan(iS /%)) 



\w + x\ m+1 (x 2 + w 2 ) 

\w + y\ m+1 {v - 2y]dy e n( a rct a n(w/y)) 

-w \V~ w\ m ~ 1 (w 2 -y 2 ) 



(40) 



As w — > then pi(x) = 0(1) if x < and 

px(x) - 2e}- v -^l^w 2 !{x) (41) 

for x > 0. To get the last asymptotic one takes into 
account that the main contribution in the integral in (]40[) 
is brought by the right small vicinity of y = 0, say (0, s), 
and apply 

exp(— aTctan(y / 6) / 5)dy ~ / cxp(—y/6 2 )dy^5 2 . 
io Jo 

(42) 

for all s > and small 8. 

The continuity of the solution of Eq. (fT6|) in the limit 
W — » is therefore guaranteed by the conditions 



C 



Ci 



Co = 2Cw 2 



(43) 



This conditions are indeed equivalent to the normaliza- 
tions conditions 



p(x)dx = 1, / q{x)dx = 



(44) 



where q = — (C + (l/4r 2 + s /t — x 2 )p). 

To show this we notice that Eq. (|T6|) can be written as 
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--(^-^ ) )P + Ox 2 ) +uq = (45) 



It follows that 



J q(y)dy = ~ ~ (472 ~- t2 ) j (p( x ) ~ p(~ x )) 

(46) 



8 



Thus, the second condition in (|44j) is equivalent to 
lim x (p(x) — p(— x)) = 



(47) 



According to Eq. (JT^J) p(x) is written as the sum of p\ 
and p 2 , whose asymptotic behavior is 



1 Pi 



and 



p+ p- 
P2(x) ^, x -> oo, _p 2 (a;) ^, a; 



(48) 



-oo (49) 



Notice that P 2 ~ = 0, P 2 + = e (."-^">/^ . Thus (f¥7|) 
becomes 



C 2 



p + 



In the limit w — > 

Pf = o(P+), P+ ~ 2^ 2 P+ 

and we get 

C 2 - 2w 2 C 



(50) 



(51) 



(52) 



which implies the continuity of the pdf. Normalization of 
p(x) in the limit w — > is equivalent to the first condition 
in 63}. 
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